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Abstract. Elastic electromagnetic nucleon form factors have long provided vital 
Q I information about the structure and composition of these most basic elements of 

' nuclear physics. The form factors are a measurable and physical manifestation of 

the nature of the nucleons' constituents and the dynamics that binds them together. 
Accurate form factor data obtained in recent years using modern experimental facilities 
has spurred a significant reevaluation of the nucleon and pictures of its structure; e.g., 
^ , the role of quark orbital angular momentum, the scale at which perturbative QCD 

(•~^ ' effects should become evident, the strangeness content, and meson-cloud effects. We 

lO ' provide a succinct survey of the experimental studies and theoretical interpretation of 

nucleon electromagnetic form factors. 
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1. Introduction 

^ I It might be said that the impact of information on nucleon electromagnetic form factors 
was first felt in 1933 when Otto Stern measured the proton's magnetic moment: 

., = (i+iJ2)5^. (1) 

The deviation from one expressed by the underlined term within the parentheses 
announced that the proton is not a point particle. Were the proton point-like, it would be 
described by Dirac's theory of relativistic fermions and hence have a magnetic moment, 
/ip = /iD = e/(2M). 

The impact continues to the present day, with modern, high-luminosity 
experimental facilities that employ large momentum transfer reactions providing 
remarkable and intriguing new information on nucleon structure [1, 2, 3]. For examples 
one need only look so far as the discrepancy between the ratio of electromagnetic proton 
form factors extracted via Rosenbluth separation and that inferred from polarisation 
transfer [4, 5, 6, 7] and the evolving picture of the contribution of s-quarks to the 
proton's electric and magnetic form factors [8, 9, 10, 11, 12, 13]. 
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1.1. Nucleon electromagnetic form factors 

In a Poincare covariant treatment an on-shell = nucleon with four-momentum 
P, mass M and polarisation s is described by a four-component spinor (column-vector) 
u{P, s) that satisfies the Dirac equation: 

(f - M)u{P, s) = 0. (2) 

In the nucleon's rest frame, s'^ = (0, s*) with s ■ s = 1. The adjoint spinor is 
u{P,s) = m(P, s)^7° and we choose to normahse such that u{P, s')u{P, s) = 2M6s's- 
(NB. Our Minkowski space metric and Dirac matrix conventions are those of [14, 15].) 
The electromagnetic current for such a nucleon is 

J,{P', s'; P, s) = u{P', s') A^(g, P) u{P, s) , (3) 

= u(P', s') (l,F,{q') + ^ ta,. F^{q^)) u{P, s) , (4) 

where P, s (P', s') are the four- momentum and polarisation of the incoming (outgoing) 
nucleon and q = P' — P is the momentum transfer. The quantity A^(g, P) in (3) 
is the fully-amputated nucleon-photon vertex, which in QCD is an eight-point Green 
function. Poincare covariance entails that the general expression for A^(g, P) involves 
twelve independent scalar functions but when projected via on-shell nucleon spinors, as 
in (3), all tensor structures reduce to the two shown in (4) with the Poincare-invariant 
scalar coefficient functions Pi and P2, which are termed, respectively, the Dirac and 
Pauli form factors. It is sometimes useful to work with the isoscalar and isovector 
combinations of these form factors 

F: = l{Ft + Fn, F: = \{Ff-Fr), (^=1,2). (5) 

Two important combinations of the Dirac and Pauli form factors are the so-called 
Sachs form factors [16]: 

Ge{Q') = F,{Q') - ^F^iQ') , Gm(Q') = Pi(Q') + F,{Q') (6) 

(— = > defines spacelike), which express the nucleon's electric and magnetic 
form factors. In the Breit frame and in the nonrelativistic limit, the three-dimensional 
Fourier transform of GEiQ"^) provides the electric-charge-density distribution within 
the nucleon, while that of GAf(Q^) gives the magnetic-current-density distribution [17]. 
Naturally, G^^{0) = 1, G^i^O) = 0, which expresses the proton and neutron electric 
charges. Moreover, Gm(0) = (Ge(0) + k) =: fi defines the proton and neutron magnetic 
moments. In this expression k = P2(0) is the anomalous magnetic moment, discovered 
for the proton by Stern: = 1.79 in (1). It is noteworthy, as we said, that for a point 
particle P2 = 0, in which case Ge = Gm- 

As reviewed, e.g., in [1, 3], it is possible to expose the contribution of individual 
quark fiavours to these form factors by considering the coupling of the Z°-boson to the 
nucleon. This is expressed via the nucleon's neutral weak current 

J^^(P', s'- P, s) = u{P', s') Aj(g, P) n(P, s) (7) 
= u{P', s') (7.Ff (g^) + ^ FiiQ') + 757.G^(Q^) + Tsf Gp(g^)) «(P, .) , (8) 
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where the new scalar functions appearing in (8) are the axial-vector, Ga-, and 
pseudoscalar, Gp, nucleon form factors, associated with the axial- vector part of the 
Z°-nucleon coupling, which herein we cannot discuss further. 

The contribution from each quark flavour to a given form factor is defined by writing 
that form factor as a sum in which each of the terms is multiplied by the relevant electric 
or weak quark charge; viz., for the proton's electric and magnetic form factors, 

GUAQ") = I g'eM') - \ \GtMiQ') + G'lM')] , (9) 



3 

" 3 

-(l-^sin^^H^) 



GtMm = {^~U^'Ow)G^^^,,{Q'' 



G^MiQ') + G'lM') , (10) 



where 9w is the weak-mixing angle: siia^ 9w{Mz) ~ 0.23. Here the form factors 
are the same in each line because, once the charges are factorised, the matrix 
elements are constructed from the same quark-level vector current whether the probe is 
electromagnetic or weak-vector. The contribution from quarks heavier than strange is 
supposed to be small. 

If one assumes charge symmetry then the (i-quark contribution to the neutron's 
form factors is the same as the u-quark contribution to the proton's; i.e., G]f^^ = G^j^j^j, 
G^M — G^EMy s-quark contribution is the same in both the proton and neutron; 

viz., G^E^^j = Giy.^ = G% j,j. Hence 

Gl^Mm = (1 -^sm'ew)GlM') - Gl,,{Q') - G%^M% (H) 

so that the s-quark contributions to the electric and magnetic form factors are 
determined once one has G^, G" and G^ . The latter is accessible via parity violating 
electron scattering from the proton, which is covered at length in [3]. 

Form factors are truly an important gauge of a hadron's structure. They are a 
measurable and physical manifestation of the nature of the hadron's constituents and 
the dynamics that binds them together. In analogy with optics, the scattering of a probe 
with three-momentum Q should resolve the structure of the hadron on a length-scale 
d ~ Hence, so long as \Q\ ^ M; namely, recoil effects are small, then one has 

a straightforward interpretation of the data in terms of a static charge- and current- 
distribution [18, 19, 20, 21, 22] and, e.g., the proton's electric and magnetic radii are 
determined via 



Ge,m{^) 



■ Ge,m{Q'^ 



(12) 

Q2=0 



dQ 

[NB. Should G(0) = 0, this normalising factor is omitted; e.g., (17) and (25).] Even 
though these aspects of the interpretation break-down for > M^, the form factors 
remain as a visible probe and exacting test of our understanding of QCD's dynamics. 

1.2. Electron scattering formalism 



Electron-nucleon scattering is typically treated in the single-photon-exchange 
approximation. As intrinsic properties of the target, the nucleon form factors are 
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independent of this approximation but their inference from experiment is not. In the 
single-photon-exchange approximation the differential cross section for elastic electron- 
nucleon scattering can be expressed by the Rosenbluth [23] formula 

^ = aM(Fl{Q^)+TFi{Q')+2T[F,{Q')+F,{Q')fi^r.\dj2)) , (13) 



Ti 



Gim+eGim , (14) 



where r = (5^/(4M^), ctm = i^^^^^j^)^ is the Mott cross-section for scattering 
from a point-like particle, uqed is the fine-structure constant, E^. is the initial electron 
energy, and E'^ and 6e are the scattered electron energy and angle, respectively. 
Reference [24] re-expressed the Rosenbluth formula in terms of Ge and Gm'-, viz., 
dc a M 
In ~ e(l + r) 

where £ = !/[! + 2(1 + r) tan^(6'e/2)] is the virtual photon polarisation parameter. 

With polarised electron beams, one also can measure the cross section asymmetry 
from a polarised target, or the polarisation transfer to an unpolarised nucleon. For a 
polarised nucleon target, the beam-target asymmetry is [25, 26]: 

Ki cos 6* GIj + K,2 sin 6* cos (f)*GEGM ni cos 6* + K2 sin 9* cos (j)*R 
^ Gl + ^sGl, ^ ^+^3 ' ^ ^ 

where ^1^2,3 are kinematic factors, 6*, 0* are the polar and azimuthal angles between 
the scattering plane and nucleon spin direction, and R = Ge/Gm- Note that in this 
case one can vary 6* to isolate the transverse component {6* = 90°), which is very 
sensitive to R, or the longitudinal component {6* = 0), which has little sensitivity to 
the form factors unless R^ ^ /^s. One can also measure the transverse and longitudinal 
polarisation transfer components which have similar sensitivity to R. 



1.3. Layout 

We have divided this article into five sections. This closes the Introduction. Section 2 
presents an experimental perspective, with an overview of the current status of nucleon 
form factor measurements. Section 3 provides a snapshot of the impact that these 
measurements are having on the phenomenology and theory of hadron physics and 
QCD. The challenges that a description of nucleon form factors poses for numerical 
simulations of lattice-regularised QCD and its successes are summarised in Section 4. 
Section 5 is an Epilogue. 



2. Experimental status 

Until very recently, most extractions of the nucleon elastic form factors came from 
unpolarised, inclusive electron scattering measurements. While this provided a 
significant body of data on G^, measurements of G^e were limited at high Q^, and 
measurements of the neutron form factors were much less precise. 

In the last decade, high intensity beams, large and efficient neutron detectors, 
and high polarisation beams and targets have led to a dramatic improvement in our 
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knowledge of the electromagnetic form factors. These new experiments, coupled with 
a better understanding of higher order electromagnetic corrections (specifically, two- 
photon exchange effects), have led to a dramatic increase in the kinematic coverage, 
precision and completeness of the form factor data. In particular, there will soon be 
measurements of all four elastic electromagnetic form factors up to ~ 3.5 GeV^, 
allowing for comparisons with nucleon structure models at both low Q^, where the pion 
cloud is believed to play an important role, and large Q^, where one is sensitive to the 
nucleon's quark core. In addition, comparisons of the proton and neutron form factors 
enable one to make model-independent statements about flavour-dependent aspects of 
the nucleon. Since contemporary simulations of lattice-regularised QCD cannot handle 
disconnected diagrams, this also provides access to the isovector form factors for which 
lattice calculations are currently feasible. 

2.1. Unpolarised elastic and quasi-elastic scattering 

Prior to the advent of polarisation transfer measurements, our knowledge of the 
proton form factors came almost entirely from inclusive elastic e-p scattering, with 
the form factors extracted using the Rosenbluth separation technique that relies on the 
simple dependence of the reduced cross-section, aR, on the virtual photon polarisation 
parameter, e [see (14)]: 



allowing one to extract Gm from at e = and Ge from the e dependence. These 
extractions have several important limitations; viz., reduced sensitivity to G^; at large 
and G^ at small (except for 6 180°), a large anti-correlation between the errors 
in the extracted values of G^; and Gm, and a significant correlation between the values 
of G^; (or Gm) extracted at different from a single experiment owing to uncertainties 
in the e-dependent corrections applied to the data. 

Despite these limitations, data from such measurements were sufficient to provide 
high precision extractions of G^ for ^ [0.1,30]GeV^ and G% from 0.01-2 GeV^ 
(see Figure 1). The uncertainties on G^ grow rapidly with owing to the reduced 
contribution of Gg to a^. Both form factors are reasonably well approximated by a 
simple dipole form, G^ = GlJup = 1/(1 + Q^/0.7lf, with in GeV^, although the 
G^j data are precise enough to show clear deviations from this dipole fit. The G^; data 
are systematically below the dipole fit on a sizeable domain near 0.1 GeV^, but the data 
at large do not exhibit deviations of the magnitude seen in G\.j. 

Similar measurements of the neutron form factors are extremely limited. Extracting 
G^ from unpolarised scattering is nearly impossible because the small value of G^ 
provides at most 5-6% of the total e-n cross-section. While this makes extraction 
on G^ easier, the absence of a free neutron target is still a significant limitation. 
Most experiments extracted G^j using inclusive quasi-elastic electron-nucleon scattering 
from deuterium [28, 29, 30, 31, 32, 33, 34]. This requires subtracting the contribution 




(16) 
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Figure 1. Ratio of and G^j to the dipolc form: l/(l+(5^/m|,) 
Tlie hollow circles and solid curve show the global analysis of the cross-section data 
from [6]. The solid squares are the results from [7], and the crosses are the high-Q^ 
measurements from [27], which were not included in [6]. 



coming from quasi-elastic e-p scattering, yielding significant uncertainties that are 
correlated with the knowledge of the proton form factors. Such extractions also require 
corrections for the nuclear effects in deuterium. One experiment extracted from a 
coincidence D{e,e',n) measurement [35]. By tagging the struck neutron, one removes 
the contamination from e-p scattering, but becomes sensitive to knowledge of the 
neutron detection efficiency and the much larger nuclear effects involved when detecting 
the neutron as well as the electron. These experiments provided measurements of 
with ~3% uncertainties on < 0.5 GeV^, ^ 5% uncertainty for 0.5-5 GeV^, and 10- 
15% uncertainty on 5 < <10 GeV^. 

For G'^, only upper limits could be set by such measurements. At very low Q^, 
neutron-electron scattering measurements have been used to extract the neutron charge 
radius [36], which is defined via 

(rl) = -Q^.GliQ') , (17) 

because G^{0) = 0, cf. (12). Prior to recent polarisation measurements, the best 
extraction of at finite came from analyses of elastic e-D scattering, where 
a model-dependent extraction of G"^ can be obtained. This extraction is extremely 
sensitive to the N-N potential used in calculating the deuterium wave function. Thus, 
while these measurements give some indication of the dependence of G^, the values 
extracted in different analyses vary by up to a factor of two [37] . 




Figure 2. Values of taken from ratio measurements on deuterium and polarised 
■^He measurements. The circles are extractions from the ratio of e-n to e-p quasielastic 
scattering and the open triangles are from measurements on polarised '^He. The right 
panel includes high data: the solid squares are the CLAS preliminary results [43], 
and the crosses [34] and asterisks [44] indicate points taken from quasielastic e-n 
scattering on light nuclei. The solid line is a fit from [18]. 



2.2. Polarised elastic scattering and ratio measurements 

The difficulties in extracting form factors from unpolarised scattering, owing to lack of 
a free neutron target and reduced sensitivity to Ge at high Q^, made it necessary to 
employ new techniques. These techniques had been known for some time, but required 
improved beams and/or detectors, which have only recently become routinely available. 

Many techniques had been proposed to make improved measurements of the nucleon 
form factors. For GJ^, measurements of D{e,e',n) could be used to exclude proton 
scattering contributions, while a comparison of D{e,e',n) with D{e,e',p) significantly 
reduces the importance of nuclear corrections [38]. Such measurements require a 
combination of high luminosity and large or efficient neutron detectors. These have 
become available in recent years and several such measurements have been performed 
[39, 40, 41, 42, 43]. These experiments yield significantly more precise measurements of 
G^, reducing the uncertainties by roughly a factor of two up to = 1 GeV^. In the 
ratio D{e, e' , n) / D{e, e' , p) , nuclear corrections largely cancel and subtraction of a large 
proton contribution is not required. 

Polarisation measurements have also provided improved extractions of G^j. For 
polarised e-n scattering, both the parallel and perpendicular asymmetry are sensitive 
only to the ratio Ge/Gm in scattering from a free nucleon. However, the technique can 
be used to extract G\[ in quasi-elastic scattering from polarised ^He [45] . In a simplified 
picture, neglecting nuclear effects and assuming unpolarised protons, the scattering is a 
combination of scattering from two unpolarised protons and one polarised neutron. The 
parallel beam-target asymmetry, 6** = in (15), from the neutron is very well known 
because it depends only on the kinematics and a small correction from (G^/G^)^. Since 
the neutron asymmetry is well known, the experiment is essentially a measurement of 




Figure 3. Values of taken from [48]-[58], compared to the Galster parametrisation 
[37]. The right panel includes projected uncertainties for future high-Q^ measurements: 
a recently completed JLab Hall A measurement (solid circles), and an approved Hall 
C measurement (hollow circle). 



the dilution of the asymmetry coming from the contribution of the unpolarised protons. 
Thus, it is essentially a measurement of the ratio of den to dep. As with the direct 
ratio measurements on deuterium, this yields minimal sensitivity to the proton form 
factors. The analysis of these experiments [46, 47] takes into account nuclear effects, 
the polarisation of the protons, and all other effects neglected in the simple picture. 
Figure 2 shows the extracted values of G^j from ratio measurements on deuterium and 
the quasielastic asymmetry in polarised '^He. 

Owing to their sensitivity to Ge/Gm, measurements of spin-dependent e-n 
scattering, via polarisation transfer or beam-target asymmetry measurements, have 
been used more extensively to extract G^. The perpendicular asymmetry, 6* = 90° 
and 0* = in (15), is sensitive to the ratio Ge/Gm [25, 26]. It provides a better way to 
measure G^ than the Rosenbluth separation technique, where the relative contribution 
of is suppressed as (G^/G^//)^- Measurements of the asymmetry have been made 
using quasielastic scattering from polarised deuterium [48, 49] and ^He [50, 51, 52, 53] 
targets, as well as measurements of recoil polarisation from an unpolarised deuterium 
target [54, 55, 56, 57]. Figure 3 shows the results of the G^ measurements described 
above. Also shown (as crosses) are the > 0.4 GeV^ results from a more recent 
extraction of via deuterium form factor measurements [58]. 

Finally, polarisation measurements have also been used to improve high-Q^ 
measurements of G^;, where the cross-section is dominated by the contribution from 
G^j. Two measurements used the cross-section asymmetry on polarised proton targets 
[59, 60], but most measurements, including all the high-Q^ extractions, have used 
recoil polarisation [4, 5, 61, 62, 63, 64]. Along with taking advantage of the increased 
sensitivity to G^; in the transverse polarisation transfer (or perpendicular asymmetry), 
some experiments have made simultaneous extractions of the longitudinal and transverse 
polarisation transfer. By taking the ratio of these components, one maintains sensitivity 




Figure 4. Polarisation measurements of jipG^/G^, showing the combined statistical 
and systematic uncertainties. The recoil polarisation data are shown as solid triangles 
[61], squares [63], stars [64], and hollow circles [4, 5, 62]. The polarised target 
measurements are shown as hollow stars [59], and crosses [60]. The solid line is a 
linear fit to the JLab data, ^pG^/G^^ = 1 - 0.13(Q2 _ o.04) [5]. 



to G^^/G\j while cancelling the uncertainties associated with the beam polarisation and 
the analysing power of the recoil polarimeter. Combining measurements with positive 
and negative helicity electrons, false asymmetries in the recoil polarimeter also cancel, 
yielding an extraction of G^^/G\j with extremely small systematic uncertainties [4]. 

Results from the polarisation transfer measurements of G^/G^j^ are shown in 
Figure 4. These measurements indicate a significant decrease in the ratio at large and 
suggest that may possess a zero at somewhat higher Q^. This differs markedly from 
conclusions based on Rosenbluth separation extractions, which indicated that G^^/G^j^ 
was roughly consistent with unity up to ~ 6 GeV^. Notably, the precise data from 
BLAST [60] have somewhat smaller uncertainties than a global analysis of \ow-Q^ cross- 
section data [20]. Hence, while in a combined analysis at low-Q^ of the cross-section 
and polarisation data the extraction of and G^j will remain limited by normalisation 
uncertainties in the cross-section measurements, the addition of the BLAST data will 
significantly reduce the anti-correlation between the extracted values of Gg and G^. 

When the discrepancy in Figure 4 was first observed, it was noted that there is 
significant scatter in the values of Gg extracted from different Rosenbluth measurements, 
especially at values above 1-2 GeV^, which is not seen in the combined analysis 
depicted in Figure 1. Owing to this scatter it was often assumed that the discrepancy 
between Rosenbluth and polarisation measurements arose from systematic uncertainties 
in the Rosenbluth extractions. However, the scatter is largely explained by the fact that 
many experiments extracted the form factors by combining new cross-section data with 
results from previous measurements. Differences in the relative normalisation of low- 
er and high-e data sets leads to significant errors in the extracted values of G^ that 
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affect tlie entire range of tlie analysis. Tlie normalisation uncertainties were often 
ignored entirely, or else estimated but treated as uncorrelated. It was demonstrated 
[65] that the various data sets were consistent when performing Rosenbluth separations 
using only data from individual experiments, or when properly taking the normalisation 
uncertainties into account. More recently, the discrepancy was demonstrated much more 
clearly by a new, high-precision Rosenbluth extraction of G^^/G\j [7]. By detecting 
the struck proton, rather than the scattered electron, this experiment had significantly 
smaller e-dependent corrections, allowing for smaller systematic uncertainties. 

Thus, the inconsistency implies either an error in the polarisation measurements, a 
common systematic error in the analysis of the cross-sections, or a fundamental flaw in 
one of the two techniques. Even if one assumes that the problem lies with the Rosenbluth 
data, owing to the greater sensitivity of the polarisation measurements, it is important 
to fully understand the discrepancy. This because the cross-section data are needed to 
separate Ge and Gm, and both the elastic cross-section and form factor measurements 
are important ingredients in other analyses. The source of the difference will determine 
its impact on other measurements [6]. Currently, the discrepancy is believed to be 
explained by two-photon exchange corrections, which are discussed in Section 3.5. 

2. 3. Future form factor measurements 

In the next few years final results should become available from the CLAS 
measurements at JLab [66], the G^ experiments in Hall A [67] and the BLAST 
experiment at MIT-Bates. Additional data on G\j were taken using a deuterium target, 
with a low momentum spectator proton tagged to give a nearly- free neutron target [68]. 
An extension of the polarisation transfer measurement of G^^/^^m — 8-5 GeV^ will 

be performed in 2007 [69], along with a series of high-precision Rosenbluth extractions of 
G^E and [70]. In the longer term, the energy upgrade at Jefferson Lab will allow high 
precision measurements of G^^l^^i G'^ to ~ 14 GeV^, as well as an extension 
of G% measurements to 8 GeV^. 

2.4- Flavour decomposition, parity violating scattering 

The combination of proton and neutron form factors can be used to study the flavour 
dependence of the charge and magnetisation distributions. Combining proton and 
neutron measurements provides sensitivity to the difference between up and down quark 
distributions because of the difference in relative weighting between the proton and 
neutron. To capitalise on this, one needs data for G^ and G% (or G\^ and G\j) that 
cover the same range and have comparable precision. With the recent and upcoming 
measurements of neutron form factors, one will have measurements of all four form 
factors at values up to ~5 GeV^. The neutron measurements have larger relative 
uncertainties, but the absolute uncertainties for the proton and neutron measurements 
are comparable. This provides an optimal case for comparison of neutron and proton 
measurements, which can be used to study the flavour dependence. 
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Adding measurements of parity violating elastic scattering allows one to fully 
separate the m, d, and s quark contributions, (7)-(ll). However, in addition to 
requiring this new data, a careful treatment of the correlations between the form factor 
measurements is required as well as calculations of the two-photon-exchange (TPE) 
corrections. Section 3.5. While the size of TPE corrections in the parity violating 
asymmetry is small [71], the impact of TPE in extracting the electromagnetic form 
factors is large enough that is must be accounted for in a combined analysis of parity- 
conserving and parity- violating elastic scattering [72], especially for >1 GeV^. 



3. Theoretical understanding 
3.1. Background 

The experimental results reviewed in Section 2 are the objective facts for which strong 
interaction theory must provide an understanding. To this we now turn. 

At large spacelike Q^, perturbative QCD amplitudes factorise. It follows that 
helicity is conserved at leading twist and hence for = —q^ > CpQCD ^ ^qcd 

p^^Qi-^ ~ constant ; (18) 

viz., the Pauli form factor is power-law suppressed with respect to the Dirac because 
it is an helicity-fiip amphtude [73]. Indeed, dimensional counting rules for QCD's hard 
amplitudes give [74, 75, 76] F^Q^) ~ l/Q^ F^^Q^) ~ l/g^ and hence 

«^>Cp^QCD ^^g^ 



Perturbative QCD and dimensional counting rules cannot predict the value of CpQCd; 
namely, the scale at which these results should become evident in experiment. That 
requires a nonperturbative method. However, the experiments which form the basis of 
Figure 4 suggest strongly that CpQCD ^ 6 GeV^. Moreover, a linear fit to the polarisation 
transfer data yields a zero in G\[Q'^^ at ^ 7.8 GeV^, while G\j{Q'^) remains 
positive definite. The possibility of a zero in G^^{Q'^) was largely overlooked before 
the polarisation transfer data became available, although thereafter many models were 
found to exhibit such behaviour. The confirmation of a zero is sought in a forthcoming 
JLab experiment [69], which will obtain polarisation transfer data out to = 8.5 GeV^. 

While a zero in G^^i^Q"^) was not generally anticipated, a ratio Rp{Q^) = 
lipG^^iQ'^) / G^j^Q"^) that falls with increasing Q"^ could early have been inferred from 
vector meson dominance fits to existing data [77]. In a later dispersion-relation fit to 
electron-nucleon scattering cross-sections [78] the ratio was explicitly calculated with an 
aim to testing (19) and found to fall with increasing on a domain 0.5 < ^^(GeV^) < 
3.0. This behaviour persists in more recent analyses of this type [79, 80, 81]. 
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3.2. Mean- field and potential models 

At its simplest, the nucleon is a nonperturbative three-body bound-state problem, 
an exact solution of which is difficult to obtain even if the interactions are known. 
Hitherto, therefore, phenomenological mean-field models have been widely employed 
to describe nucleon structure; e.g., soliton models [82, 83, 84] and constituent-quark 
models [85, 86, 87]. Such models are most naturally applied to processes involving small 
momentum transfer {Q^ < M^, M is the nucleon mass) but, as commonly formulated, 
their applicability may be extended to processes involving moderately larger momentum 
transfer by working in the Breit frame [88]. An alternative is to define an equivalent, 
Galilean invariant Hamiltonian and reinterpret that as the Poincare invariant mass 
operator for a quantum mechanical theory [89, 90]. This is likely a truer extension of a 
model to large-Q^. 

In the context of soliton models, a zero in G^(Q^) at ~ 10 GeV^ was evident 
in a study of nucleon electromagnetic form factors based on a topological model that 
incorporated coupling to vector mesons and relativistic recoil corrections [91]. However, 
owing to the uncertainties noted above in treating relativistic recoil corrections, it was 
argued therein that the predictive power of the model was poor at large-Q^. 

Appropriate to the study of electromagnetic form factors at > M^, [92] employed 
a light-front constituent- quark model with a wave-function parametrisation introduced 
in [93, 94]. In an earlier application, this model was found to produce a proton electric 
form factor with a zero at ~ 5.7 GeV^ [95]. Although the predicted position of the 
zero is contradicted by existing data, it is a feature of this type of model that a zero in 
G^{Q^) is easily obtained. 

The data in Figure 4 can be re-expressed as a measure of KpF2 (Q^) / Ff {Q'^) . This 
ratio is approximately constant on 2 < (^^(GeV^) < 6. Such behaviour has been argued 
to indicate the presence of substantial quark orbital angular momentum in the proton 
[19, 92, 96, 97, 98]. It is a feature of [92] that VQ^ Fi (Q^) / F[ (Q^) ^constant for 
2 ^ Q^(GeV^) < 20. However, such behaviour is not uniformly found. Moreover, it 
has been argued [99] that the dependence of form factors on y/Q^ instead of is an 
artefact of the construction employed in [93, 94]. The lack of analyticity in can be 
related to a violation of crossing symmetry. Crossing symmetry is a necessary property 
of quantum field theory. Its violation in the model of refs. [92, 93, 94, 95, 100] indicates 
that while the model is relativistic and may provide an efficacious tool for computing 
low-Q^ properties of hadrons, it is not applicable at large momentum transfer because 
it is inconsistent with essential features of quantum field theory. 

A comparative analysis of relativistic constituent-quark model calculations of 
baryon form factors is provided in [101]. Model results depend in general on the 
representation of the baryon mass operator, which is expressed in [101] via algebraic 
parametrisations of the ground-state wave function. They also depend on the form of 
the current-operators, which in [101] are generated by the dynamics from single-quark 
currents that are Poincare covariant under a kinematic subgroup: i.e., instant-form. 
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point-form [102] or front-form [92, 95, 99, 103]. (NB. The point-form construction of 
spectator-current operators is not fully constrained by Poincare invariance alone [104]. 
One can appeal to time-reversal invariance to further limit the composition, but the 
result still contains a residual arbitrariness which increases with [105].) Notably, 
the studies in [102, 103] are predictive and falsifiable because they are based on a mass 
operator from which wave-functions are actually calculated. 

Instant-form kinematics emphasises features that are closest to those of 
nonrelativistic quark models, with an interpretation of the wave-function that stresses 
covariance under three-dimensional rotations and translations. In models based on 
light-front kinematics the relevant Lorentz boosts and translations are all kinematic, 
whereas in the point-form approach there are no kinematical translations and also no 
kinematic interpretation of the wave function as a representation of spatial structure. In 
a relativistic approach the addition of spins to form the total angular momentum of the 
composite system always requires momentum-dependent rotations of constituent spins 
[106]. These are the Wigner rotations (related to Melosh rotations in the front-from), 
which generate a purely kinematical contribution to constituent-quark orbital angular 
momentum and guarantee that the bound-state wave-function is an eigenfunction of J 
and Jz in its rest frame. 

With an aim of providing a broad perspective, [101] illustrated that, independent 
of the form of kinematics used in constructing the current operator, a simple algebraic 
representation of the ground state orbital wave function can yield a fair description 
of individual elastic form factors. The proton's electric form factor was found most 
sensitive to the form of kinematics and the corresponding spectator current. The 
best-fit parameters that characterise the wave-function vary markedly from one form 
of kinematics to another. These parameters determine the range and shape of the 
orbital wave-function. In order to fit the data, instant- and front-form kinematics 
demand a spatially extended wave function, whereas a compact wave-function is required 
in point-form kinematics. To quantify, using one measure of the rms matter radius 
of the constituent-quark wave-function, the instant- and front-form wave-functions 
correspond to tq ^ 0.6 fm, whereas the point-form wave-function is characterised by 
ro ~ 0.2 fm. In connection with Figure 4, and as illustrated in Figure 5, instant- 
form dynamics as implemented in [101] gives a ratio that exhibits modest suppression, 
matching most closely the form factor data inferred via Rosenbluth separation in the one- 
photon exchange approximation. The point-form calculation drops more rapidly with 
increasing Q^; and the front-form result drops most rapidly, lying below but tracking the 
polarisation transfer data. It is thus plain that accurate form factor data can be used 
to test constituent-quark model mass operators. Equally, however, such mass operators 
are model-specific and do not possess a veracious connection with QCD dynamics. 
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Figure 5. MpG^/G^j vs. in relativistic quantum meclianics: Solid-curve - instant 
form; dotted-curve - point form; dashed-curve - front form. For context, the open 
symbols indicate the ratio inferred from JLab polarisation transfer data and the filled 
triangles are from Rosenbluth separation. (Figure adapted from [101].) 



3.3. Three-body problem in quantum field theory 

In quantum field theory a meson (quark- ant iquark bound-state) appears as a pole 
in a four-point quark-antiquark Green function (see, e.g., [107], and [108] for 
phenomenological applications). The residue of that pole is proportional to the meson's 
Bethe-Salpeter amplitude, which is determined by an homogeneous Bethe-Salpeter 
equation. It is plain by analogy that a nucleon (three-quark bound-state) must appear 
as a pole in a six-point quark Green function. (These facts underpin the extraction 
of masses and form factors through numerical simulations of lattice-regularised QCD, 
Section 4.) In this case the residue of the pole is proportional to the nucleon's Faddeev 
amplitude, which is obtained from a Poincare covariant Faddeev equation that adds-up 
all the possible quantum field theoretical exchanges and interactions that can take place 
between the three dressed- quarks that constitute the nucleon. 

While this is true in principle, the tractable treatment of the Faddeev equation in 
quantum field theory requires a truncation. One such scheme is founded [109, 110] on 
the observation that an interaction which describes colour-singlet mesons also generates 
quark-quark (diquark) correlations in the colour-3 (antitriplet) channel [111]. Naturally, 
diquarks are confined; namely, they are not asymptotic states in the strong interaction 
spectrum [112, 113, 114, 115, 116]. Reference [117] reported a rudimentary study of this 
Faddeev equation and subsequently more sophisticated analyses have appeared; e.g., 
[118, 119, 120, 121]. It has become apparent that the dominant correlations for ground 
state octet and decuplet baryons are scalar and axial-vector diquarks. This may be 
understood on the grounds that: the associated mass-scales are smaller than the masses 
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of these baryons [122, 123], with models typically giving masses (in GeV) 

"^Mo+ = ^•'^^ ~ ^■^'^ ' = = ^{dd)^+ = 0.95 - 1.02 ; (20) 

the electromagnetic size of these correlations is less than that of the proton [124] 

r[„,]^^ ^0.7fm, (21) 

from which one may estimate ~ 0.8 fm based on the p-meson/vr-meson radius- 

ratio [125, 126]; and the positive parity of the correlations matches that of the baryons. 
Both scalar and axial- vector diquarks provide attraction in the Faddeev equation; e.g., 
a scalar diquark alone provides for a bound octet baryon and including axial-vector 
correlations reduces that baryon's mass. 

For some it may be noteworthy that the possibility of diquark correlation has been 
studied in numerical simulations of lattice-regularised quenched-QCD. Reference [127] 
obtains so-called weakly bound scalar and axial-vector diquarks whose masses are in 
accord with (20). This is confirmed in subsequent analyses that are not uniformly 
restricted to the quenched-truncation; e.g., [128, 129].^ In addition, a lattice estimate of 
diquark size [128] is consistent with (21). 

The truncation of the Faddeev equation's kernel is completed by specifying that the 
quarks within the baryon are dressed, with two of the three dressed- quarks correlated 
always as a colour-3 diquark. Binding is then effected by the iterated exchange of roles 
between the bystander and diquark-participant quarks. This ensures that the Faddeev 
amplitude exhibits the correct symmetry properties under fermion interchange. A Ward- 
Takahashi-identity-preserving electromagnetic current for the baryon thus constituted 
is subsequently derived [130]. It depends on the electromagnetic properties of the 
axial-vector diquark correlation: its magnetic and quadrupole moments, and 

respectively; and the strength of electromagnetically induced O''" ^ 1"^ diquark 
transitions, nr- 

Thus does one arrive at an analogue in quantum field theory of the mass- and 
current-operators necessary in the quantum mechanical treatments described above. 
A merit of the Poincare covariant Faddeev equation is that a modern understanding 
of the structure of dressed- quarks and -gluons is straightforwardly incorporated; viz., 
effects owing to and arising from the strong momentum dependence of these propagators 
are realised and exhibited. This momentum- dependence, which explains, e.g., the 
connection between constituent- and current-quark masses, was predicted by Dyson- 
Schwinger equation studies and has been confirmed in lattice simulations. A synopsis 
can be found in Section 5.1 of [131]. 

Figure 6 depicts the calculated [133, 134] ratio of the proton's Sachs form factors; 
viz., /ipG^(Q^)/G^((5^). It is important to understand the behaviour of the data cf. 

f It is curious that some refer to the O"*" correlation as the good diquark and the l"*" correlation as 
the bad diquark. We emphasise that the axial-vector diquark enables dynamical correlations within a 
baryon that are essential. A realistic description is forfeited if the axial-vector diquark is omitted. 




Figure 6. fJ'pG^/G^^ vs. Q^. was calculated using the point-particle values: 

= 2 and xi+ ~ li K-r = 2. Variations in the axial- vector diquark parameters 
used to evaluate have little effect on the results. The width of the band reflects the 
variation in G^j with axial-vector diquark parameters: the upper border is obtained 
with Hi+ = 3, Xi+ = 1 f-iid KT ~ 2, while the lower has = 1. Data: squares-l^]; 
diamonds-[5]; triangles-[7] and circles-[132\. (Figure adapted from [133] and prepared 
with the assistance of A. H611.) 



the calculation at small Q^. In the neighbourhood of = 0, 

GliQ^) ^i_Q^ \(r ?- (r^"?] (22) 

where Vp and are, respectively, the electric and magnetic radii. Experimentally, 
Tp ^ and this explains why the data varies by less than 10% on < < 0.6 GeV^. 
The calculated curve was obtained ignoring the contribution from pion loops, which 
interfere constructively with those from the axial-vector diquark correlations. Without 
such chiral corrections, Vp > and hence the calculated ratio falls immediately with 
increasing Q^. Incorporating pion loops one readily finds Vp ^ [134]. It is thus 
apparent that the small behaviour of this ratio is materially affected by the proton's 
pion cloud. Moreover, such contributions may actually play a role to < 2 GeV^. 

These features of pseudoscalar meson contributions to form factors are evident, 
e.g., in [100, 135, 136, 137] and we judge that a pointwise accurate description of the 
individual proton and neutron form factors at small is impossible without a careful 
treatment of meson cloud effects. This is particularly true of the neutron's charge 
form factor, since G^{Q^ = 0) = and the form factor is never large, so its evolution 
is sensitive to delicate cancellations between the quark-core and meson cloud. Note, 
however, that merely perceiving a small-Q^ deviation between the form factors and a 
global dipole fit is not an unambiguous signal of pion cloud effects because a priori 
there is no reason to expect the form factors to be accurately described by a single 
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dipole parametrisation valid uniformly on a large Q^-domain. 

Pseudoscalar mesons are not pointlike and hence pion cloud contributions to form 
factors diminish in magnitude with increasing . It follows that the evolution of 
IXpG^^iQ"^) / G^j^Q"^) on Q"^ > 2GeV^ is primarily determined by the quark core of the 
proton. This is evident in Figure 6, which illustrates that for >2 GeV^, /ip G^^/G\.j is 
sensitive to the parameters defining the axial-vector-diquark-photon vertex. The ratio 
passes through zero at ^ 6.5 GeV^; namely, at the point for which (^^(Q^) = 0. 
In this approach the existence of the zero is robust but its location depends on the 
model's parameters. The behaviour of UpG^^/G^j owes itself primarily to spin-isospin 
correlations in the nucleon's Faddeev amplitude. The forthcoming JLab experiment [69] 
will test these predictions. 

Naturally, the ratio HnG^{Q'^)/G'^{Q'^) is also of experimental and theoretical 
importance. Notably, in the neighbourhood of = 0, 

GUQ') ~ & ^ ^ 

where r„ is the neutron's charge radius. The Faddeev approach shows (23) to be a good 
approximation for r^Q^ < 1 [138] and extant data [56] are consistent with this. It is thus 
evident that, as for the proton, this ratio's small behaviour is materially affected by 
the neutron's pion cloud. Reference [138] predicts that the ratio will continue to increase 
steadily until — 8 GeV^. That will be examined in future experiments. Section 2.3. 

In common with relativistic constituent-quark models, the Faddeev equation 
analysis described above makes an assumption about the dynamical content of QCD. 
The assumptions can be wrong but they are testable. Moreover, in this case the 
Schwinger functions at which one arrives can in principle; viz., at some future time, 
be calculated via numerical simulations of lattice-regularised QCD. 

The predictions of perturbative QCD were revisited in an analysis [139] that 
considers effects arising from both the proton's leading- and subleading-twist light-cone 
wave functions, the latter of which represents quarks with one unit of orbital angular 
momentum, with the result 



[lnQ7A2]2+^ 



constant, > A , (24) 



where /3 = 11 — |iVj, with Nf the number of quark flavours of mass ^ Q^, and A 
is a mass-scale that corresponds to an upper-bound on the domain of nonperturbative 
(soft) momenta. This is naturally just a refined version of (18). Equation (24) is not 
predictive unless the value of A is known a priori. However, A cannot be computed 
in perturbation theory. Notwithstanding this the empirical observation was made [139] 
that on the domain 2 < (^^(GeV^) < 6 the Q^-dependence of the polarisation transfer 
data is approximately described by (24) with 0.2 < A(GeV) < 0.4. However, that 
is merely accidental: A ~ 0.3 GeV corresponds to a length scale ta ~ Ifm and it is 
not credible that perturbative QCD is valid at ranges greater than the proton's radius. 
In fact, one can argue [134] that a judicious estimate of the least-upper-bound on the 




Figure 7. Proton Pauli/Dirac form factor ratio, calculated using a Poincare covariant 
Faddeev equation model [134] and expressed via (24) with A = 0.94 GeV. The band is 
as described in Figure 6 except that here the upper border is obtained with = 1, 
Xi+ = 1 and Kq- = 2, and the lower with = 3. The data are: squares-^i] and 
diamonds-[5]. (Figure adapted from [133] and prepared with the assistance of A. H611.) 



domain of soft momenta is A = Mj^; viz., the nucleon mass, in which case the Dirac 
and Pauh form factors obey (24) for > 5 GeV^ as illustrated in Figure 7. 

The result in Figure 7 is not significantly influenced by details of the diquarks' 
electromagnetic properties. Instead, the behaviour is primarily governed by correlations 
expressed in the proton's Faddeev amplitude and, in particular, by the amount of 
intrinsic quark orbital angular momentum [98]. The nature of the kernel in the Faddeev 
equation (or, analogously, a mass operator) specifies just how much quark orbital angular 
momentum is present in a baryon's rest frame. 

It is noteworthy that orbital angular momentum is not a Poincare invariant. 
However, if absent in a particular frame, it will inevitably appear in another frame 
related via a Poincare transformation. (Therefore is Wigner rotation necessary in 
constituent-quark models.) Nonzero quark orbital angular momentum is a necessary 
outcome of a Poincare covariant description, which is why the covariant Faddeev 
amplitude is a matrix- valued function with a rich structure that, in a baryons' rest frame, 
corresponds to a relativistic wave function with s-, p- and even d-wave components. This 
is well illustrated by Figure 6 in [119], which explicitly depicts these components of the 
Faddeev wave function in the nucleon's rest frame. A crude estimate based on their 
magnitudes indicates that the probability for a u-quark to carry the proton's spin is 
P„t ~ 80 %, with P„| ~ 5 %, ~ 5 % and Pdi ~ 10 %. Hence, by this reckoning 
~ 30% of the proton's rest-frame spin is located in dressed-quark angular momentum. 

Elastic electromagnetic form factors can be expressed as one- dimensional integrals 
of valence-quark generalised parton distributions (GPDs), which are nonforward matrix 



Nucleon electromagnetic form factors 



19 



0.15 
0.1 
0.05 



-0.05 
-0.1 
-0.15 

-1.5 -1 -0.5 0.5 1 1.5 



Figure 8. The contours display the 68 and 95% confidence intervals from an analysis 
[143] of G^" and G^J at = 0.1 GeV^, together with estimates of these quantities 
from model calculations [144, 145, 146, 147, 148, 149, 150] and numerical simulations 
of lattice-regularised QCD [151, 152, 153, 154]. Recent JLab measurements [13] not 
included in [143] are overlaid for comparison: diagonal band - HAPPEx-H; horizontal 
hand - HAPPEx-^Ie. (Figure composed with assistance of R. D. Young.) 




elements of light-front operators. It follows that accurate experimental measurements 
and theoretical calculations of electromagnetic form factors can be used to place much- 
needed stringent constraints on parametrisations of GPDs. With such parametrisations 
in hand one can, e.g., estimate the contribution of quark spin and orbital angular 
momentum to the light-front nucleon spin [96, 140, 141, 142]. 

3.4- Strangeness in the proton 

As we described briefly in Sections 1.1 and 2.4, the s-quark contribution to the proton's 
form factors is accessible via parity violating electron scattering if one has accurately 
determined G^^j^j{Q^), G^^^(Q^). Naturally, since the nucleon has no net strangeness, 
G^^{0) = 0. However, there is no such simple constraint on either the sign or magnitude 
of = G^(0). In analogy with (17), a strangeness charge-radius can be defined via 

i<s) = -6 ^-^j^ ■ (25) 

In Figure 8 we provide a snapshot of the current status of experiment and theory 
for the strange form factors of the proton. One model estimate lies within the 
95% confidence limit [149]. It is inferred from a dispersion-relation fit to nucleon 
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Figure 9. New version of a global analysis of the L-T (Rosenbluth) data on fipCs/ Gm 
[6]. Open circles - ratio obtained from cross-sections before two-photon correction; 
closed circles - the same ratio obtained after correcting the cross-sections for two- 
photon exchange [159]. Dashed line: a fit to the high-Q^ polarisation transfer data [6]. 



electromagnetic form factors and yields /i^ = 0.003, r^^ = 0.002 fm^. The lattice-QCD 
estimates are described in Section 4. It has been argued [143] that at present the world's 
data are consistent with the strange form factors of the proton being zero. Forthcoming 
experimental results from JLab and Mainz will contribute more to the picture that is 
developing of the form factors accessible through parity violating electron scattering. 

3. 5. Two-photon exchange: Rosenbluth and polarisation transfer 

It is apparent from our discussion that theory views the polarisation transfer data as the 
truest measure of G^/G^j. There is a mounting body of evidence that the discrepancy 
between this data and that obtained via Rosenbluth separation arises from the effects of 
two-photon exchange (TPE) contributions to the cross-section, which are only partially 
accounted for in the standard treatment of radiative corrections [155]. For example, it 
was demonstrated [156] that TPE contributions could explain the discrepancy without 
spoiling the linear ^-dependence of the reduced cross section. Moreover, improved 
evaluations of the effect of the exchange of an additional soft photon [157, 158] had 
already indicated that the effects were larger than previously believed. 

Figure 9 illustrates the impact that TPE corrections can have. Calculations 
including hard and soft photons have recently been performed in both hadronic 
[159, 160, 161, 162] and partonic [163] models. However, at this time there is no complete 
calculation, valid at all relevant kinematics. Hadronic calculations generally account 
only for an unexcited proton in the hadronic intermediate state. They are thus limited 
to relatively low where the contributions of the excited intermediate states should 
be small. Contributions beyond the elastic, including the A [160] and higher mass 
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states are under investigation. The partonic calculations are only valid at large s, t and 
values, and rely on a model of the proton's GPDs. One can also make empirical 
estimates of the TPE amplitudes [156, 164] but this requires assumptions about the 
angular dependence of the amplitudes and is limited by the precision of the Rosenbluth 
data when using the discrepancy as a measure of the size of the TPE effects. 

These calculations and empirical estimates of the TPE amplitudes make several 
common predictions. First, all can at least partially resolve the discrepancy between the 
Rosenbluth and polarisation transfer results, with TPE mainly changing as extracted 
from Rosenbluth separations. They all point to relatively weak Q^- dependence, meaning 
that the corrections could also impact upon precision measurements at low Q^. Finally, 
they indicate that the effects are largest at large scattering angles, corresponding to 
e — s> 0, and small for 5 — 1, which implies that there will also be an impact on the 
extraction of G\,[ from the extrapolation of aRto e = 0, (16). 

The idea that TPE corrections could be large enough to explain the discrepancy 
seemed originally to contradict limits set by previous measurements designed to 
test the Born approximation. First, the linear e-dependence of the reduced cross- 
section is consistent with one-photon exchange and second, several comparisons of 
e'^-p and e~-p (and fi^-p) scattering showed no indication of TPE effects (see, e.g., 
[155, 165, 166, 167, 168, 169]). The interference between Born and TPE amplitudes 
changes sign when the sign of the charge of one of the scattered particles is flipped, and 
these measurements indicated that TPE effects were small, on the order of 1%, while 
the discrepancy implies missing corrections of up to 6% [6]. Finally, measurements 
of the normal-polarisation in polarisation transfer experiments [170, 171, 172] and the 
asymmetry. An, in beam-target asymmetry measurements [173, 174], are consistent 
with zero, in accordance with the Born approximation expectation. 

It has since been found that each of these earlier measurements had insufficient 
sensitivity or kinematic coverage to observe TPE effects of the scale predicted by 
contemporary calculations. Owing to low luminosity of the secondary positron and 
muon beams, almost all the comparisons were at very low Q^, or at small scattering 
angles, corresponding to e > 0.7. The limited data at smaller e, on the domain where 
calculations indicate that TPE contributions are largest, provide evidence for a TPE 
correction of (4.9 ± 1.4)% [175]. Similarly, while a global analysis was able to set tight 
limits on the size of nonlinearities in the e-dependence of the reduced cross section 
[176], these limits are consistent with the calculated nonlinearities [159, 163]. New 
positron measurements will be performed in the near future to better measure the 
effect of TPE at large angles and moderate values [177, 178]. Further high- precision 
Rosenbluth separation measurements [70], using the proton detection technique of [7], 
will provide the sensitivity needed to observe the calculated nonlinearities. It will also 
provide clean measurements of TPE effects on the cross-section at large values, where 
the contribution of becomes negligible. 

At present the form factor discrepancy provides indirect evidence for TPE, while the 
existing positron measurements at large scattering angle provide only a 3a indication. 
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While future measurements will significantly improve the situation, there is one 
observable that has provided clear evidence for the presence of TPE corrections in elastic 
e-p scattering. The asymmetry in the scattering of transversely polarised electrons from 
unpolarised protons, Aj^, is zero in the Born approximation but can be non-zero owing to 
TPE contributions. Measurements have shown significant asymmetries in measurements 
at MIT-Bates [179] and Mainz [180]. Both Ax_ and A^r are sensitive to the imaginary part 
of the TPE amplitude, while the form factors depend on the real part of the amplitude. 
Hence these measurements do not provide direct quantitative information that can be 
applied to the form factor measurements. However, they do provide a clean signal for 
the presence of TPE and can be used to test calculations of the TPE amplitudes. 

There is an ongoing experimental and theoretical program aimed at fully 
understanding TPE corrections. Whereas much of the focus has been on the discrepancy 
in G^; at large Q^, TPE corrections are also important at low and in extracting G\j. 
Since TPE corrections vary slowly with their effects will also be important for high- 
precision measurements at low . While the fractional TPE correction to is large, 
it is typically comparable to or only slightly larger than the experimental uncertainties. 
The correction to G\.^ is much smaller, typically at the few percent level, but this is 
larger than the quoted experimental uncertainties. Hence the correction to G\.^ can 
potentially have a larger impact in constraining calculations of the form factors. 

The quantitative study of the effects of TPE has only just begun. The impact 
of these corrections at low has been evaluated in extractions of the proton charge 
radius [21] and the s-quark contribution to the nucleon form factors [72]. However, these 
form factors are important input for many analyses. The effects of TPE corrections 
must therefore be more quantitatively understood so that we may have the reliable 
input needed to study; e.g., hyperfine splitting in hydrogen [181], L-T separations in 
quasielastic scattering [182], extractions of the axial form factor from neutrino scattering 
[183], and VCS/DVCS measurements [184] that require precise knowledge of the Bethe- 
Heitler process in order to extract the Compton scattering amplitudes. It is not always 
sufficient to provide corrected form factors, as many of these processes will also be 
modified by TPE and the corrections will not be identical for all observables. 

4. Lattice QCD 

Section 2 explains that nucleon form factors have been studied experimentally to very 
high precision, while Section 3 exhibits that this precision data is available on the 
nonperturbative domain of QCD whereupon models and truncations of QCD have been 
extensively applied. On this domain the numerical simulation of lattice-regularised 
QCD, when used in conjunction with chiral effective field theory, is becoming an 
important addition to the array of tools available in modern hadron physics. Moreover, 
it is widely hoped that at some future time this approach may yield uniquely reliable 
predictions for hadron observables from QCD. 
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4--1- Overview 

The lattice study of hadron electromagnetic form factors has a long history. The 
first calculations of the pion form factor were performed over 20 years ago, initially 
with SU{2) colour [185, 186, 187] and subsequently 5f/(3) colour [188, 189]. They 
were promptly followed by calculations of the proton's electric form factor [190]. Soon 
thereafter, magnetic moments and electric charge radii were extracted from the tt, p and 
N electromagnetic form factors [191]. Despite the limited computing resources of the 
late '80s and early '90s, it proved possible to perform calculations of the electromagnetic 
properties, including magnetic moments and charge radii, of the entire baryon octet 
[192] and decuplet [193], with results that could be compared with experiment. The 
first attempt at examining the Q^-dependence of nucleon electromagnetic form factors 
was reported in [194], where the authors found positive values for G^. 

There has recently been renewed interest in calculating electromagnetic form factors 
on the lattice. For example, the QCDSF collaboration performed a quenched-QCD 
analysis of these form factors at momentum transfers G (0.45, 1.95) GeV^ [195]. 
(NB. The domain is fixed and restricted by the lattice regularisation. According to 
Section 1.1, such momentum transfers resolve length-scales d ~ 0.14 — 0.3 fm. The 
proton's charge radius is = 0.85 fm.) The study was performed at three quark masses 
and three lattice spacings, allowing the chiral and continuum limits to be investigated. 
Upon comparison with experiment, the authors observed the lattice dipole masses to 
be too large [see the discussion of (26)] and attributed this to quenching. However, the 
limited spatial resolution might also bear some of the responsibility. 

This analysis was extended in [196] to include an extraction of magnetic moments 
and charge radii. That allowed for a comparison with chiral effective field theory 
(ChEFT). The authors derived an extrapolating function from ChEFT that was 
then compared with the lattice results, which were calculated at heavy quark masses 
(m^ > 500 MeV). Their extrapolations of the proton and neutron anomalous magnetic 
moments confirmed predictions [197, 198, 199, 200, 201, 202] that substantial curvature 
is required between the lowest-m^r lattice-data point and the chiral limit in order to 
obtain agreement between the lattice results and experiment. This substantial curvature 
was also predicted in [203] , where quenched lattice results were extrapolated to the chiral 
limit using finite range regularisation [204]. 

Hitherto the most extensive study of was carried out in quenched-QCD [205] 
with Q"^ G (0.3, 1.0) GeV^. (This corresponds to a resolution of length-scales d ~ 0.20 
- 0.36 fm. The scale associated with the neutron radius is 0.58 fm.) The authors 
observed a positive form factor, to which a Galster parametrisation [37] was fitted 
and therefrom a value inferred for the neutron charge radius at each of the simulated 
quark masses. Various methods of extrapolating to the chiral limit were subsequently 
considered. It was argued that the radius deduced could be reconciled with experiment 
by standard phenomenology and lowest- or next-to-lowest-order contributions from 
chiral perturbation theory. A result obtained more directly is lacking. 
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The renewed focus on form factor calculations may be attributed to the challenge 
of new precision data, and improvements in both algorithms and machine speed, which 
are enabling simulations to be performed with dynamical quarks, at quark masses 
corresponding to ^ 300 GeV, on lattices with spacing a ~ 0.1 fm and spatial extent 
L ~ 2.5 fm. This brings the promise of more realistic simulation results. It is nonetheless 
sobering that the computational demands are estimated to increase as 



nQ') = 7rTT^i-TV2 ' (26) 



4-. 2. Magnetic Moments and Form Factor Radii 

As explained in association with (12), radii are extracted from a form factor's slope 
at = 0, and magnetic moments are obtained from Gm{Q^ = 0) [see the discussion 
associated with (6)]. However, lattice kinematics entails that = is not directly 
accessible in the simulations. Hence an extrapolation from 7^ is required in order 
to infer values of these static properties and one commonly fits the form factor results 
using a dipole: 

with -F(O) the fitted normalisation and tud the fitted dipole mass. 

A large statistics investigation of the electromagnetic properties of the octet baryons 
in quenched-QCD has recently been performed [206]. Magnetic moments, and electric 
and magnetic radii were extracted from the form factors for each individual quark- 
fiavour in order to test the environmental sensitivity of the quark contributions to 
these quantities. Simulations were performed with pion masses as low as 300 MeV 
in order to search for evidence of the chiral nonanalytic behavior predicted by quenched 
chiral perturbation theory. Of particular interest was an observed environmental isospin 
dependence of the strange quark distributions in A° and S'^. It was found that when the 
environmental quarks are in an isospin-0 state (A baryon) the strange quark distribution 
is broader than when the environmental quarks are in an isospin-1 state (E baryon). 

An up-to-date study of isovector nucleon electromagnetic form factors, (5), at 
momentum transfers G (0.2,2.5) GeV^ is reported in [207]. This kinematic domain 
probes length-scales in the range (0.12, 0.4) fm. The calculations were performed in 
quenched-QCD and in Nf = 2-QCD. At the pion masses accessible in this study; i.e., 
^0.4 GeV, the effects of unquenching were perceived to be small - an observation 
consistent with the analysis of [203]. In comparison with experiment, the lattice results 
for the isovector form factors in [207] again lie uniformly above the data. For the 
ratio Rp described in Section 3.1, the lattice results are constant out to = 2.5 GeV^ 
whereas the polarisation transfer data show suppression (see Figure 4). With a chiral 
extrapolation used to estimate the isovector magnetic moment in the chiral limit, 
agreement with the experimental value was obtained. On the other hand, consistent 
with earlier calculations, the charge radius is constant on the domain of quark masses 
employed and hence disagrees with experiment. It is possible that the mismatch with 
experiment owes to current-quark masses that are too far from reality and/or finite 



Nucleon electromagnetic form factors 



25 



0.8 
0.6 
0.4 
0.2 





experiment 
400 MeV 
550 MeV 
700 MeV 



• • • 



0.5 



1.5 2 

[GeV]' 



2.5 



3.5 



0.6 
0.5 
0.4 
0.3 
0.2 
0.1 






A 


Nf=2 [208] 




A 


Nf=0 [196] 







JNf=U yiy) 1 J 




♦ 


Nf=2 [207] 








f 




** * & » - 









0.2 



0.4 



0.6 

[GeV^] 



0.8 



1.2 



Figure 10. Left panel: F\ form factor for three different pion masses compared with 
experimental data [208]. Right panel: Isovector F\ (Dirac) form factor radius as a 
function of from both quenched-QCD [196, 207] and iV/ = 2-QCD [207, 208]. 



lattice spacing effects. Reference [207] noted that if dynamical fermions and small quark 
masses are needed, then a realistic study will require large computer resources. 

Aspects of the behaviour just described are illustrated in Figure 10. The left panel 
depicts a typical example of the proton's Dirac form factor calculated at three different 
pion masses [208]. The results sit high cf. experiment, with a slight trend towards the 
experimental points as m^r is decreased. This translates into small form factor radii 
that increase with decreasing pion mass, as seen in the right panel, wherein we display 
quenched-QCD [196, 206, 207] and iV/ = 2-QCD [207, 208] results for the isovector Fi 
form factor radius, plotted as a function of m^. The figure indicates a gradual evolution 
toward the chiral limit. However, ChEFT predicts that both the Fx and F2 radii should 
increase dramatically in the neighbourhood of the chiral limit [196, 203]. In modern 
studies there are indications of incipient "chiral curvature" for < 400 MeV. Recent 
results from the LHPC Collaboration [209] also exhibit the features described here. 

The left panel of Figure 11 depicts the latest lattice results for the isovector 
magnetic moment (inferred by fitting (26) on the accessible kinematic domain and listed 
in units of lattice nuclear magnetons, iin = e/2M''^", where M''^" is the nucleon mass 
measured in the lattice simulation) as a function of m^, obtained using both quenched- 
QCD [196, 206, 207] and Nf = 2-QCD [207, 210]. The experimental value is indicated 
by a star at the physical pion mass. Consistent with [203], there appears to be little 
difference between the results obtained in quenched and unquenched simulations. The 
evolution with appears gentle. However, it is plain that a linear extrapolation would 
miss the experimental point by many standard deviations. This, too, is compatible with 
chiral perturbation theory, which predicts a rapid increase in /i^r as —>■ [196, 203]. 
(There is a hint of this chiral curvature at the smaller m.,^ values in Figure 11.) 
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QCD [207, 210]. Right panel: F2/F1 form factor ratio at « 550 MeV and three 
different lattice spacings [208, 210] plotted according to (24) with A = 0.2 GeV. 



4-3. Investigating the dependence and reaching for large 

As reviewed in Section 3.1, counting rules for QCD's hard amplitudes predict Fi(Q^) ~ 
l/Q^ and F2{Q'^) ~ It is currently difficult, however, to obtain lattice results with 

high enough precision over a large enough range of values to distinguish between 
dipole and quadrupole behaviour. Furthermore, it is extremely likely that the present 
limitations on lattice kinematics prevent a determination of CpQCD and an exploration 
of the domain Q"^ > CpQCD ^ ^qcd which such behaviour should be apparent. 

Some attempts have nevertheless been made in connection with the ratio 
F2{Q'^) / FilQ"^) , whose perturbative scaling is described by (24). The right panel of 
Figure 11 displays values for the ratio in (24) calculated with A = 0.2 GeV from results 
at three lattice spacings with approximately the same pion mass [208, 210]: the ratio 
is roughly constant on the domain of explored. This notwithstanding, as discussed 
in connection with (24) and Figure 7, such a low value of A is not credible as the 
least-upper-bound on the domain of soft momenta in QCD. 

On the qualitative side, the results in Fig. 11, together with those in [207], indicate 
that lattice spacing, quark mass and quenching effects on F2(Q^)/Fi(Q^) appear small 
in comparison with statistical errors. Quantitatively, though, the lattice values are 
higher than the corresponding experimental data. Thus, while contemporary lattice 
simulations are able to reproduce qualitative features of the experimental data, smaller 
pion masses, at least, are needed before a quantitative description can become a reality. 

Herein we have focused attention on the Q^-evolution of the ratio fipG^^{Q'^) / G\j{Q'^) 
because of the intriguing possibility that it will pass through zero at some > 5 GeV^. 
In Section 4.2 we highlighted Nf = 2-QCD studies [207] that explore nucleon form fac- 
tors for G (0.2, 2.5) GeV^ and remarked that the ratio reported therein is constant 
(or perhaps grows slightly) with increasing in marked contrast to JLab's polarisa- 
tion transfer data. Figure 4. There is a hint in the lattice results that with the lightest 
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accessible current-quark mass the ratio dips below one at ~ 1.5 GeV^. This provides 
encouragement for lattice simulations planned at more realistic quark masses. 

While simulations with lighter quark masses are beginning to become possible, 
the problem remains that a typical lattice calculation is restricted to a small domain 
of relatively low momenta. This has inspired the LHPC Collaboration to assess the 
computation cost for a calculation of the electromagnetic form factors out to ~ 
6 GeV^ [211]- They focused on the Breit frame, [p' = —p = 2TTn/L), because previous 
studies indicated that therein the form factors have smaller statistical uncertainties than 
with other momentum combinations at the same Q^. Their analysis reveals that the 
relative error in F^IQ'^), (5), at fixed pion mass increases as n^. Since their point with 
= 4 (Q^ = 4.15 GeV^) has a relative error of 62%, then in order to achieve a point 
at = 8 {Q"^ ~ 6 GeV^) with a relative error of 30%, they would have to increase the 
statistical accuracy by at least a factor of 50. Furthermore, to compound the difficulty, 
it was observed that the relative error in the isovector Dirac form factor increased with 
approximately the fourth power of l/m^r. It is evident therefore that an essentially new 
technique must be found before numerical simulations of lattice-regularised QCD are in 
a position to calculate the form factors at large with realistic quark masses. 

4.4- Strangeness content of the nucleon 

The determination of the strange quark content of the nucleon offers a unique 
opportunity to obtain information on the role of hidden flavour in the structure of 
the nucleon (see Sections 1.1, 2.4 and 3.4). 

Direct lattice-QCD calculations of the strangeness content are computationally 
demanding in the extreme and results have so far proved to be inconclusive [151, 152, 
212]. While there is optimism that increases in computing power may enable the next 
generation of lattice calculations to obtain accurate results, this may actually require 
an investigation into new lattice techniques. One possibility is the background fleld 
method [213, 214, 215, 216], where a weak signal may be enhanced by coupling a 
strong electromagnetic fleld to the vacuum strange quarks. This technique has recently 
been employed successfully to calculate physical quantities, such as magnetic moments 
[217], and electric [218] and magnetic [219] polarisabilities. Alternatively, a method 
of evaluating the all-to-all propagator, developed by the Dublin group [220], offers 
signiflcantly improved precision over traditional stochastic estimators, and it would be 
interesting to see this employed in a strangeness form factor calculation. 

Meanwhile, one must continue to rely on more indirect methods for an extraction of 
the strangeness form factors. Constraints of charge symmetry were combined with chiral 
extrapolation techniques, based on flnite-range-regularisation, and low-mass quenched- 
QCD lattice simulations of the individual quark contributions to the charge radii and 
magnetic moments of the nucleon octet, to obtain precise estimates of the proton's 
strange electric charge radius [154] and magnetic moment [153]: 

(r2)P = +0.001 ± 0.004 ± 0.002 fm^ , /x, = -0.046 ± 0.019 /x^v • (27) 
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Together, these resuhs considerably constrain the role of hidden flavour in the structure 
of the nucleon. Moreover, Figure 8 shows that they agree extremely well with a recent 
analysis [143] of the world's complete data set on parity violating electron scattering. 
They are also consistent with the latest measurements from Jefferson Lab [13], which 
are indicated by the coloured bands in the figure. 

5. Epilogue 

The world's hadron physics facilities are providing data of unprecedented accuracy, 
from both nuclear and hadronic targets. Herein we have focused on nucleon elastic 
electromagnetic form factors but, before closing, a short digression on the pion is 
worthwhile. In that case, too, the impact on our understanding of the basic features of 
QCD is enormous. 

The pion is notionally a two-body bound-state and hence the simplest composite 
system described by QCD. However, the pion is also QCD's Goldstone mode and its 
properties are intimately connected with confinement and dynamical chiral symmetry 
breaking (DCSB). This dichotomy can only be reconciled, and a veracious explanation 
of pion properties thereby achieved, through a treatment using the full machinery of 
quantum field theory, and the approach employed in this must guarantee that the 
axial- vector Ward-Takahashi identity is accurately realised [221]. DCSB is expressed 
through this identity, and there are strong indications that DCSB is a necessary 
consequence of confinement. Thus, information on the pion probes into the deepest 
part of QCD. New data [222] and the reanalysis of old data [223] are confronting theory, 
e.g. [224, 225]. Through this, reaching to higher momentum transfers promises to 
identify those elements fundamental to a precise understanding. Moreover, there is 
room for some optimism that the JLab upgrade will enable data to be acquired on a 
domain [226] in which there are hints of the transition to truly perturbative behaviour. 

Returning to nucleon elastic electromagnetic form factors, the next few years 
will see new neutron data and an extension of proton electric form factor data to 
= 8.5 GeV^, while the middle of the next decade should see precision neutron and 
proton data to ~ 14GeV^. This is part of a much larger programme that will 
yield an enormous body of concrete information about the spectrum of hadrons and the 
interactions between them. An accurate understanding of this data will draw a map of 
the distribution of mass and spin within the nucleon; lay out the connection between 
the current-quark and the constituent-quark; and should enable us to become certain 
of the domain on which perturbative QCD becomes predictive. 
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